Enhanced dynamical entanglement transfer with multiple qubits 
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We present two strategies to enhance the dynamical entanglement transfer from continuous vari- 
able (CV) to finite dimensional systems by employing multiple qubits. First, we consider the entan- 
glement transfer to a composite finite dimensional system of many qubits simultaneously interacting 
with a bipartite CV field. We show that, considering realistic conditions in the generation of CV 
entanglement, a small number of qubits resonantly coupled to the CV system is sufficient for an 
almost complete dynamical transfer of the entanglement. Our analysis also sheds further light on 
the transition between microscopic and macroscopic behaviours of composite finite dimensional sys- 
tems coupled to bosonic fields (like atomic clouds interacting with light). Furthermore, we present 
a protocol based on sequential interactions of the CV system with some ancillary qubit systems and 
on subsequent measurements, allowing to probabilistically convert CV entanglement into 'almost 
perfect' Bell pairs of two qubits. Our proposals are suited for realizations in various experimental 
settings, ranging from cavity-QED to cavity-integrated superconducting devices. 

PACS numbers: 03.67.Mn,42.50.Dv,42.50.Pq 



I. INTRODUCTION 



Due to the infinite dimension of their Hilbert space, 
continuous variable (CV) systems allow, in principle, to 
store an infinite amount of entanglement. Even in re- 
alistic settings, the entanglement of two-mode squeezed 
states generated in parametric amplification processes 
with present technology can be much more than one ebit, 
thereby largely exceeding the maximum amount of entan- 
glement between a pair of qubits. Nonetheless, discrete 
variable systems - and, most notably, qubits - are natu- 
rally privileged for the implementation of many quantum 
information tasks, with the prominent example of quan- 
tum computation. In this respect, it is of interest to 
investigate the efficiency at which entanglement can be 
transferred from continuous to discrete variable systems 
(by realistic coherent interactions) and to envisage strate- 
gies to improve such efficiency. This is relevant in view of 
the relative ease with which highly entangled CV states 
can be currently generated With efficient trans- 

fer procedures, such entanglement could be distributed 
to separate discrete variable systems and employed for 
general quantum information processing purposes. 

Internal levels of atoms coherently interacting with 
bosonic light fields are natural candidates as discrete 
variable receivers of CV entanglement. In fact, the en- 
tanglement transfer between two radiation modes and 
a pair of two-level atoms throiigh coherent interactions 
has already been investigated H, ^|^. Also, the state 
transfer between macroscopic atomic clouds and light has 
been theoretically considered recently and important 
steps have been made towards its experimental demon- 
stration , relying on measurements and coherent inter- 
actions. In the limit of a macroscopic number of polarised 
atoms, the atomic component of such systems behaves 
as a CV system, so that the entanglement transfer re- 



lates, actually, two CV systems of different nature. For a 
resonant Jaynes-Cummings coupling ^J, the interaction 
reduces to the action of a beam splitter between the two 
systems allowing, in principle, perfect state (and thus 
entanglement) transfer. However, if the receiving system 
is really a discrete variable one, like a small number of 
atoms {i.e. a "microscopic" cloud) would be, then the 
transfer is no longer perfect. Some questions of funda- 
mental and practical interest arise in this instance: How 
many atoms are actually needed in order to realise an 'es- 
sentially perfect' transfer, thus exhausting the resources 
of the CV system? How does the transition between 
microscopic and macroscopic entanglement transfers be- 
have? 

The purpose of this paper is twofold. On the one 
hand, we aim to study the entanglement transfer from 
two light modes to a pair of atomic ensembles made up 
of a small number of two-level atoms (or, equivalently, to 
systems with Hilbert spaces of small dimension d > 2). 
We will show that, for two-mode squeezed states realis- 
tically achievable in the labs, all the entanglement can 
be extracted by a few atoms, in an essentially micro- 
scopic regime. Moreover, we will shed further light on the 
transition from the finite dimensional regime to the CV 
behaviour displayed by macroscopic ensembles, explicitly 
elucidating the algebraic reasons lying behind this transi- 
tion and the physical conditions allowing to treat macro- 
scopic ensembles as CV systems. On the other hand, we 
propose a method to increase the entanglement transfer 
to a pair of two-level atoms (qubits), by letting further 
pairs interact successively with the same entangled light 
field and then by postselecting the local measurements 
on such pairs. We will show that, following this route, 
the entanglement transferred to the first pair of qubits in- 
creases considerably with respect to the strategy pursued 
in Refs. Remarkably, we demonstrate that such an 
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'extraction' procedure, aiming at achieving a Bell state 
of first two atoms, can be made "arbitrarily perfect" by 
repeated iterations of the probabilistic protocol. 

Both the schemes we suggest in this work are natu- 
rally implemented in the context of cavity-quantum elec- 
trodynamics (cavity-QED). In this scenario, the qubits 
would be embodied by two-level atoms interacting with 
the field modes of two optical cavities externally fed by 
the two-mode CV system. The feeding process can be ac- 
curately described by the effective mixing of the external 
CV state with vacuum at a beam splitter (rem-esented by 
the cavity lossy mirror), as detailed in Refs. [^^|^. The 
trade-off between a cavity which can be externally fed 
and a sufficiently long coherence time of any intra-cavity 
photon can be accomplished by exploiting the non-linear 
effects of the intra-cavity atomic medium as described in 
Ref. ^3 (see also Ref. 0). Concerning the simultaneous 
presence of a multi-qubit system in a cavity, we recall 
that an efhcient scheme able to infer the actual number 
of atoms inside an optical resonator has been recently re- 
alised, allowing for the exact preparation of up to 3 atoms 
inside a single cavity (with the envisaged possibility of ex- 
tending this number to ~ 10), over a time three orders of 
magnitude shorter than the effective trapping time 
Because of the generality of our protocols, other physical 
systems may be considered in order to implement them. 
For instance, the postselction scheme can be suitably re- 
alised by considering systems of superconducting qubits 
integrated in quasi-unidimensional cavity-structures , 
as it will be more explicitly described in Section [llll 

This paper is organized as follows. In Sectional after 
having introduced the system at hand and the concept of 
entanglement transfer, we analyze in detail the extension 
of this protocol to multiple qubits simultaneously inter- 
acting with the CV modes. In Section lllll we address 
the case of sequential interactions of indipendent pairs of 
qubits with the CV modes, followed by proper postselec- 
tion events (occurring from the second pair on) . Finally, 
in Section [IVI we summarize our results. 



II. MULTIPLE QUBITS INTERACTING 
SIMULTANEOUSLY 



Throughout the paper, we will consider two ensem- 
bles of N two-level systems at two separate sites (we will 
indifferently refer to them as 'atoms' or 'qubits') inter- 
acting with a CV system made up of two modes of the 
radiation field (the distributor) ^ previously prepared in 
an entangled two-mode squeezed state j^^) reading 
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(1) 



where |n,n) stands for the tensor product of two num- 
ber states \n) in each mode, while r is the two-mode 
squeezing parameter. The entanglement of this state in- 
creases with increasing r, with a logarithmic negativity 



Ej^ given by Ej^f = 2r/ln(2) (see the following for the 
definition of Ej^). For simplicity, we consider ensem- 
bles of two-level atoms, each interacting with the field 
through a Jaynes-Cummings Hamiltonian with the same 
coupling strength ("Tavis-Cummings" model J^J). The 
total interaction Hamiltonian of such a system is given, 
taking into account both sites (labelled by 1 and 2), by 
Hint = Hi+ H2, with 
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where akj for k — 1,2 and j — 1, . . . ,N is the ladder 
operator referring to atom j in site k (2ct^ = (^x + i<^y hi 
term of the Pauli matrices ax and ay), while ak stands for 
the annihilation operator of the radiation mode in site k. 
Note that the calculations can be easily generalized to the 
case of asymmetric interaction which, however, would not 
significantly affect our main results. Furthermore, we will 
restrict to the case of a pure state of the CV system. A 
mixed state can be easily taken into account, and would 
simply result in the shrinking of the peak values of the 
transferred entanglement with no substantial changes 
in the qualitative behaviour of the transfer. 

We also assume that all the atoms are initialised in 
the ground state at the beginning of the interaction i.e., 
adopting the spin terminology, that they are perfectly 
"polarised". Note that this assumption, besides being 
practically satisfiable by the application of a strong po- 
larising field, is also reasonable, as it has already been 
shown that, for = 1, such an initial atomic con- 
figuration allows for an optimal entanglement transfer 
3]. Because of this hypothesis, the global state of the 
atoms (each of which can be seen as a pseudo-spin 1/2) 
will stay within the symmetric (A -I- 1) -dimensional sub- 
space with total spin j = A/2 during the evolution. 
The dynamics at each site thus reduces to the inter- 
action of one mode of the field with a spin j initially 
prepared in the ground state |j, m — — j) (m being the 
value of the projection of the total angular momentum 
vector along the quantization axis), with raising oper- 
ator cr^. The action of a'^ on an atomic state is de- 
termined by the Clebsch-Gordan coefficients according 
to = = V U -m){j + m+ 1) \j, m + 1)^, for 

m — — J, . . . , j (k — 1,2 referring to the site). For conve- 
nience, let us relabel the atomic states {\j,rn)f,} at each 
site by their number n of 'excitations' over the ground 
state: |n)j. = \j, n ~ j)^. for n = 0, . . . , A. 

The dynamics at each site can be analytically treated 
for A < 4 14]. For instance, let us recall the time evo- 
lution operator for A = 1: 
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where the matrix acts on the ordered atomic two- 
dimensional Hilbert space {|1) , |0)}fe (whereas a and 
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act on the infinite-dimensional Hilbert space of the field), 
T = gt is the rescaled interaction time and h = 1. We 
have numerically integrated the dynamics for larger num- 
ber of atoms (up to iV = 20), truncating the infinite 
dimensional Hilbert space of the field modes at num- 
ber states depending on the initial squeezing parameter 
r (obviously, as r increases higher excitation states get 
more significantly populated). For the sake of clarity, 
let us explicitly describe the dynamics of the system un- 
der the action of the evolution operator exp(— «ifi„ji) — 
exp{—iHit) (g) exp{—iH2t) = Ui®U2- The density ma- 
trix q{t) describing the state of the system at the rescaled 
time T reads 

q{t) = C/i ® U2{\^r) K |0,0))((0,0| K (e.|)C/t ® Ul , (4) 

where (8) denotes the tensor product with respect to the 
bipartition in the two sites, H denotes the tensor prod- 
uct with respect to the bipartition into the atomic and 
radiation subsystems and the state |0,0) = |0)j^ |0)2 
stands for the ground state of the global atomic system. 
We are now interested in the reduced atomic state Qat{T), 
obtained by partially tracing the global state £i(t) over 
the field variables: Qat{T) = Tr/gi(T), Tr/ denoting the 
trace over the Hilbert space of the field. 

Finally, the 'logarithmic negativity' £V of the state 
Qat{T) can be worked out to estimate the rate of en- 
tanglement transferred from the CV two-mode squeezed 
state to the finite dimensional Hilbert space of the atoms. 
Let us recall that such a quantity is strictly related to 
the positivity of the partial transpose (PPT) criterion 
for the separability of quantum states [l5j. The partial 
transposed density matrix g is obtained from any given 
bipartite quantum state g by transposing the variables of 
only one of the two subsystems. The PPT criterion then 
simply reads g>Q. Indeed, such a criterion is necessary 
for the separability of any quantum state regardless of 
the dimension of the system's Hilbert space but is also 
sufficient only for 2 x 2 and 3x2 Hilbert spaces (consid- 
ering only finite-dimensional cases). The 'negativity' M 
(first envisaged in Rcfs. later thoroughly discussed 
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0.70 
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1.10 


2.07 


0.65 




10 


1.35 


2.77 


0.71 



TABLE I: Optimal values of the two-mode squeezing param- 
eter, maximal logarithmic negativity achievable (through res- 
onant interaction with a CV state with the optimal two-mode 
squeezing) and 'efficiency' ef f of the optimal transfer [quanti- 
fied by the ratio between the transferred logarithmic negativ- 
ity and the logarithmic negativity of the initial distributor's 
CV state: eff = ln(2)/(2r)] for the first values of TV. 
Note the monotonic increase of the ratio eff. 
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FIG. 1: (color online). Entanglement transferred from a two- 
mode squeezed state with r ~ 0.86 (with initial Em — 2.48) 
for various number A'^ of receiving two-level atoms at each 
site.. The red curve refers to A'^ = 15 (solid line), the blue 
curve to A^ = 7 (dashed line), the cyan curve to A' = 2 (dash- 
dotted line) and the green curve to A = 1 (dotted line). All 
quantities plotted are dimensionless. 



in Refs. |17|). which can be easily determined from the 
knowledge of the density matrix (upon diagonalisation of 
g), is simply defined as the absolute value of the sum of 
the negative eigenvalues of g and directly quantifies the 
violation of the PPT criterion. The logarithmic negativ- 
ity is then just given by £V = log2(2A/' -I- 1) (so that 
Em' — 1 for Bell states of two qubits). From an opera- 
tional point of view, the quantity Ef^ is an upper bound 
to the distillable entanglement and is related to the en- 
tanglement cost under PPT preserving operations [T^ . 
Let us mention that, since the PPT criterion is no longer 
sufficient for Hilbert spaces of dimension 3x3 and higher, 
such a quantity is null for some entangled states in higher 
dimensions. Still, Ej^ has been recently proven to be an 
entanglement monotone irrespective of the dimension of 
the Hilbert space under study Moreover, it satis- 

fies the other crucial requirement of being null for non 
entangled states. Thus, it can be employed as a proper 
quantification of the entanglement for high-dimensional 
Hilbert spaces. Clearly, increasing the dimension of the 
Hilbert space allows for "more space" being available for 
entanglement: in fact, the maximal value of Em" is given 
by log2((i) for a bipartite (d x d)— dimensional Hilbert 
space. 

Let us now move on to analyse the entanglement trans- 
fer to pairs of "microscopic" ensembles of N atoms. In- 
terestingly, as already found in the case = 1 0, the 
maximal entanglement which can be dynamically trans- 
ferred to a finite dimensional system is not monotonically 
increasing with the entanglement of the initial CV state 
(i.e. with r). Indeed, the value of r allowing for the op- 
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(a) (b) (c) (d) 

FIG. 2: (color online). Comparison between the entanglement transfer to pair of ensembles of A'^ atoms (solid lines) and to 
pairs bosonic modes with interaction strength rescaled by \/N (dotted lines), for N=2 (a), for N=7 (b) , for N=ll (c), for 
N=15 (d). All quantities plotted are dimensionless. 



timal transfer (after maximisation over time), which we 
will denote by ropt{N), increases with increasing N but 
is always finite. The behaviour of the function ropt{N), 
whose values for small TV are reported in Table HI is not 
trivial, resulting from a compromise between the amount 
of entanglement initially present in the system (monoton- 
ically increasing with r) and a population of the excitated 
states of the field modes favouring a better state transfer 
(too large r's imply large populations in higher number 
states and, as a consequence, a mismatch of the 'effective' 
dimensions of atoms' and light's Hilbert spaces). Notice 
that 

lim ropt{N) = oo , 

A''— *oo 

as in the infinite dimensional (CV) limit a beam-splitter 
allows for perfect state transfer and, in principle, an 
infinitely entangled initial state (obtained in the limit 
r — > oo) can be also perfectly transferred to the atoms. 

Also, one has ropt(l) ~ 0.86, for which the initial CV 
state |Cropt(i)) has a logarithmic negativity £V — 2.48. 
For such a two-mode squeezing parameter, which is in 
the range of current experimental techniques 0, 01 1 we 
have investigated the possibility of realising 'exhaustive' 
entanglement transfer from CV to discrete variable sys- 
tems by increasing the number of atoms N. Fig. ^ 
shows the behavior of the transfer for different values 
of TV. The efficiency and velocity of the transfer increase 
with increasing N, the value of 15 allowing for a very 
good transfer to the 16-dimensional atomic space (with 
Ejif ~ 2.38 ~ 0. 96 Ej^). Let us mention that the trun- 
cation of the CV state to the state |15,15) has a log- 
arithmic negativity Ej^ ~ 2.47, thus containing almost 
all the entanglement of the CV state. In general, for 
the squeezing parameters currently achievable, a small 
number of atoms (TV ~ 10 — 100) is always capable to 
retrieve a substantial part of the entanglement contained 
in the initial two-mode squeezed state. For instance, the 
most entangled two-mode state reported in Ref. has 
a logarithmic negativity Ej\j- ~ 1.60, corresponding to 
r ~ 0.55 [23|. For such a value of r, a number of atoms 

= 10 allows for the transfer of a logarithmic negativity 



£V — 1.55, whereas for = 20 one reaches the value 
£V ~ 1.58, thus almost exhausting the resources of the 
initial CV state. Indeed, values of N of this order allows 
for an almost perfect extraction of the entanglement even 
in instances still far from the current experimental pos- 
sibilities. To be more specific, let us consider the most 
entangled state considered in the first of Ref. 2], with 
a logarithmic negativity E_\f ~ 4.53 corresponding to 
r ~ 1.57. Now, the logarithmic negativity of this CV 
state is mostly contained into the first 40 levels of the 
field's Hilbert space (the truncation of the CV state to 
the state |40, 40) has E_\f ~ 4.49). Even though a direct 
numerical treatment of this instance would be cumber- 
some, one can provide a very good estimate of the loga- 
rithmic negativity which can be transferred to 40 atoms 
by assuming that the same ratio between the maximal 
transferable logarithmic negativity and the initial avail- 
able one as in the case r = 0.86 and = 15 applies 
here. Note that this a very conservative assumption, as 
the efficiency of the transfer generally increases with in- 
creasing A^ since, as we will appreciate in detail shortly, 
the finite dimensional system better mimic a quantum 
harmonic oscillator (see also Tabled). Finally, we infer 
that a logarithmic negativity E_^j- > 4.31 could be trans- 
ferred from the CV to the finite dimensional system, by 
letting interact iV ~ 40 qubits simultaneously. 

As we have seen, while being a great improvement with 
respect to the transfer to single atoms, the "microscopic" 
transfer is never completely exhaustive, even though the 
receiving Hilbert space would be able to almost perfectly 
contain the initial CV state. We now show that the rea- 
son for this fact lies in the dynamics of the system, and 
not in the mismatch between the different Hilbert spaces. 
To this aim, let us consider in more detail the features of 
the transition of the atomic system between the micro- 
scopic and the macroscopic regime which, as well known, 
can mimic under certain conditions a CV system inter- 
acting through a beam-splitter. Firstly, let us remark 
that, if the ladder operators of the atomic ensembles 
were behaving like actual bosonic ladder operators b\, 
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i.e. according to 

bi\n), = V^^\n+l), (5) 

(with the additional proviso bl |iV)^. — 0), the state 
(and thus entanglement) transfer from light to the atoms 
would be essentially perfect for low enough squeezing pa- 
rameters (as the only requirement on the atomic subsys- 
tem would consist in having enough dimensions to con- 
tain the initial field's state up to its effective truncation), 
with perfect (not 'convoluted') oscillations resulting from 
a beam-splitter-like dynamics. On the other hand, the 
action of the actual operators is described by 

a+ |n), = ^{N~n){n+1) \n + 1), (6) 

for n = 0, . . . , A^. Let us now focus on the subspaces 
Efe (for k=l,2) of the atomic Hilbert spaces, spanned by 
basis- vectors |ri) j, such that N ^ n (corresponding to the 
macroscopic, "highly polarised" limit). Denoting by 
and the restrictions of cr^ and bk to the subspace Sfe 
one has, as apparent from the previous two equations, 
cr^^ ~ \/Nbl., which properly accounts for the bosonic- 
like behaviour of the polarised atomic system. Notice 
that this relationship also implies that, on the 'highly 
polarised' subspace, the SU(2) algebra approximates the 
Heisenberg algebra of the canonical commutation rela- 
tions, because the restriction of the operator can be 
regarded as a (large) constant: a^^ ~ N/2. Explicitly, 
[cr^fciO'Sfc] — N ~ ^''^Sfc' should expect. Let us 

note that the \/]V factor relating the ladder operators 
to their bosonic approximations accurately explains the 
aforementioned speed-up in the transfer with increasing 
N. To explore in detail the transition to the macroscopic 
regime, we have numerically compared the (real) entan- 
glement transfer obtained with the atomic operators 
to the ideal transfer obtained by replacing them with the 
operators . As shown in Fig. [21 the mimicking is still 
not very accurate even for N — 15, thus explaining the 
good, but not quite perfect entanglement transfer in this 
instance (even though all the needed capability would be 
available in the receiving Hilbert space). Moreover, in 
all the plots, it is apparent that the matching with the 
regular bosonic oscillations is lost with increasing inter- 
action time: this is due to the fact that higher atomic 
levels (with n not negligeable with respect to N) get ex- 
cited and the atoms lose their original polarization. In 
fact, the bosonic approximation commonly employed in 
treating the interaction of light with atomic ensembles 
is not only restricted to the macroscopic and polarised 
case, but also to short interaction times. 

As a final remark about the entanglement transfer to 
small atomic ensembles, let us mention the possibility, al- 
ready addressed in Ref . , of probabilistically improv- 
ing the performance of the transfer by measurements on 
the CV system and subsequent postselection. In partic- 
ular, we have here considered the measurement of the 
parity of the field's photon number. The results of this 
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FIG. 3: (color online). Entanglement transfer for A'" = 1 (a) 
and N — 2 (b) after postselection of the finite dimensional 
state for an even (solid line) and odd [dotted line, always null 
in case (a)] CV system's photon number. The measurements 
are taken after the CV and qubit systems have interacted for 
a rescaled time r. The non postselected, 'purely dynamical' 
entanglement transfer is plotted as a reference (dashed line). 
All quantities plotted are dimensionless. 



study for the cases A^ = 1 and N — 2 are reported in 
Fig. 13 In general, the outcome "even" always increases 
the entanglement of the finite dimensional state. Con- 
versely, the outcome "odd" always completely spoils the 
entanglement for A^ = 1 (corresponding to the tranfer 
to two two-level systems) as, because of the structure of 
the state resulting from the dynamical evolution, it ef- 
fectively projects such a state onto a separable subspace. 
Still, notice that for A^ > 1 the outcome "odd" may ac- 
tually enhance the entanglement, depending on the time 
at which the measurement is performed. 



III. MULTIPLE ATOMS INTERACTING 
SUCCESSIVELY 

As we have seen in the previous section, simultaneous 
interactions with many atoms allow for a very efficient 
entanglement extraction from infinite to finite dimen- 
sional systems. Despite this remarkable fact, the main 
interest of such entangling schemes lies in the possibility 
of achieving highly entangled states of pairs of two-level 
systems, which are privileged and archetypical for com- 
putational and communication protocols. In the present 
section we present a strategy (based on the measurement 
of ancillary qubit subsystems and on the subsequent post- 
selection) to enhance the entanglement extraction from 
the CV system to two qubits. In general, such schemes re- 
sult in the purification of the interesting two-qubit state 
and thus in a sort of entanglement 'purification', as is 
the case for the postselection of detection events of the 
field modes (see the previous section and Ref. |23)- We 
mention that a scheme employing ancillary qubits and 
measurements was also proposed in Ref. 5], addressing 
the realization of quantum repeaters through atom-light 
interactions. 

Here, instead of measuring the state of the CV sub- 
system (which might present some practical difficulties), 
we take in exam the measurement of the states of an- 
cillary qubit-pairs. The scheme of principle of the pro- 
tocol we consider here is sketched in Fig. 0] In order 
to elucidate the basics of the protocol, we consider just 
two pairs of qubits, the first one given by qubits 1 and 

2 and the second by qubits 3 and 4. The first pair in- 
teracts with the entanglement distributor for a time ti, 
after which the system is described by the density matrix 
p(ri) determined by Eq. (0J for = 1. Then, the sec- 
ond pair of qubits interacts with the CV system for an 
interaction time T2. Each qubit interacts with the local 
light mode according to the Jaynes-Cummings Hamil- 
tonian of Eq. jSJ). After such a dynamical evolution, 
qubits 3 and 4 are individually measured in the orthog- 
onal basis {\ipa) , \fpa)'^} with {ipa) = a\0) + (3 |1). Here, 

= 1 — |ap and \a\ € [0, 1] (for simplicity, we assume 
the same measurement is performed on the two qubit 
ancillary systems ^2^). We consider all the qubits ini- 
tialized in the ground state |0) and the entangled state 
l^r) of Eq. as the initial state of the field, constituting 
the entanglement resource. Let us recall, once more, that 
this choice for the initialization of the qubits, other than 
being natural, has been proven to optimize the entan- 
glement transfer process ^0,113. In what follows, we 
set the value of the two-mode squeezing parameter r to 
0.86 which, as already remarked, allows for the maximal 
transfer in the "passive" scheme with A = 1. We de- 
note by Us and U4 the evolution operators that describe 
the dynamics of qubits 3 and 4 interacting with the light 
modes [each of them being given by Eq. Q for r = T2]. 
Let us also suppose that the measurement projects qubits 

3 and 4 onto the two-qubit state \cd) = \tpa)(^\''Pa) ■ After 
such a measurement, qubits 1 and 2 are left in the state 



6 




FIG. 4: Scheme of the protocol based on sequential passage 
of a multi-qubit system, shown for just two qubit pairs. Two 
qubits, labelled 1 and 2, interact with the entanglement dis- 
tributor, represented by a two-mode squeezed state (gener- 
ated by the non-degenerate-parametric amplifier (NDPA)). 
Then, a second pair of qubits, 3 and 4, are sent to interact 
with the entangler and are subsequently measured (for in- 
stance, by means of two channeltrons) onto the qubit-basis 
{j^a) , IV'a)^} (defined in the body of the paper). 

gcd,i2 reading 

Tr34j[|cd) {cd\ (Us (g> Ui)g{n)iul (g> uj)] 
TT[\cd){cd\{U3<»U4)g{n){ul^ul)] ' 
with probability 

p = Tr [\cd) {cd\ {Us ® Ui)Q{Ti){ul ® [/])]. (8) 

Here Tr34j stands for the partial trace on the Hilbert 
space of qubits 3 and 4 and on the field while Tr is the 
usual global trace. 

Before proceeding with the quantitative description 
of the performances of this protocol, it is worth re- 
marking that the addressed situation can be realized 
in different physical setups such as cavity- QED sys- 
tems, individually-trapped ions interacting with optical 
cavities .24.1 and cavity-integrated superconducting sys- 
tems p. I12I l23j . In the first case, the two pairs of 
qubits could cross an optical cavity along two orthog- 
onal directions on the same plane (see Fig. ^25j. In 
the second and third specific qubit-pair can be 

set in resonance with a cavity-field mode by tuning an 
external magnetic fiux which, at the same time, sets the 
qubits of the second pair in a dispersive (far off-resonant) 
regime. The technology needed for this kind of oper- 
ations already exists, requiring a magnetic field with a 
spatial gradient. The magnetic field globally addresses 
the two qubits of two different pairs which are simul- 
taneously present inside a single cavity, but differently 
Zeeman-shifts them, setting one qubit in resonance and 
the other off-resonance (in a way which extends and gen- 
eralizes the Mintert-Wunderlich proposal |2^). However, 
while the strong-coupling regime seems to be remotely 
feasible in the cavity-trapped ion setup, the main ad- 
vantage which characterizes the superconducting imple- 
mentation is the natural achievement of these working 
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FIG. 5: Entanglement against the rescaled interaction time 
Ti — T for the case \a\'^ = (with both the qubits 3 and 4 
found in |1)) and an initial squeezing parameter r = 0.86. 
From (a) to (d) we show the behavior of the entanglement 
function for T2 = 0, T2 — n, T2 — 3n/2 and T2 — 2n. All 
quantities plotted are dimensionless. 



conditions. This is due to the fact that the electro- 
magnetic fields, in a system of superconducting charge 
qubits ll^l, couple directly to the excess charge of the 
qubit rather than to its intrinsic dipole moment. This 
gives rise to effective dipole moments which are orders of 
magnitude larger than the already large ones character- 
izing alkali Rydberg atoms interacting with cavity-field 
modes 0, Moreover, the very high cavity quality 

factors achievable in this quasi-unidimensional configura- 
tion 0,1131 and a judicious choice of the superconducting 
qubits working point (the charge-degeneracy point, as de- 
scribed in Refs. [3. 112.12^ 1. which strongly quenches the 
non-Mar kovian 1// noise |^ ^ , make our proposal 
robust against decoherence. Thus, a suitable modifica- 
tion of already realizable fully integrated setups seems to 
be very promisig for the implementation of the present 
proposal. 

As for the scheme of principle of our proposal, note 
that an alternative approach based on the simultaneous 
interaction of each light field with N qubits (as in the 
previous section) and on the subsequent measurements 
of — 1 qubit systems to localize the entanglement in 
the interesting pair might have been considered. Actu- 
ally, even allowing for generic qubit measurements, such 
a strategy cannot beat the direct dynamical transfer to 
two qubits (obtained for A = 1): the "dispersion" of the 
entanglement among the qubit subsystems cannot be re- 
versed by local measurements and classical communica- 
tion. The projection of qubit-pairs which sequentially in- 
teract with the distributor, as is the case for the scheme 
presented here, is thus a key to improve the 'passive', 
purely dynamical transfer. 



0.08 
0.06 
0.04 
0.02 




FIG. 6: Positive part of the difference between the logaritmic 
negativity for ti — t and T2 = 3tt/2 and the peak of logarith- 
mic negativity for ri = 3n/2 and T2 = (for which Ej^f ~ 0.9) 
against the rescaled interaction time r. We have considered 
r = 0.86 and \af = as in the previous figure. All quantities 
plotted are dimensionless. 



In Fig. |5l we show the behavior of between qubits 
1 and 2 as a function of the rescaled interaction time 
Ti = T. The interaction time T2 (relative to the cou- 
pling between qubits 3 and 4 and the CV system) is used 
as a parameter. We have considered a — {i.e. both 
the ancillary qubits are found in their excited state). As 
is apparent, the analytical structure of the transferred 
entanglement changes drastically in presence of the an- 
cillary interaction and postselection which is, not surpris- 
ingly, not always advantageous for all interaction times 
as it can even completely spoil the entanglement between 
the qubits. However, large peaks of entanglement, with 
values beating the maximum obtained with the passive 
strategy (corresponding to E_\f ~ 0.9, see Tabled, are 
frequently present. In particular, as detailed in Fig.|Sl for 
the interaction times ti = 5.65 and T2 — 3tt/2, the con- 
sidered postselection process results in £V — 0.975 with 
probability p ~ 0.3, thus allowing for the distillation of a 
quasi- maximally entangled state. 

Before addressing the refining of such almost perfect 
Bell pairs extraction procedure, some further remarks are 
in order. The procedure is qualitatively effective regard- 
less of the choice of a. Any value of a allows to beat the 
maximum of transferred entanglement achievable with 
the passive strategy (in Fig.[7|we report the analysis for 
a — 0.95; in particular, panel (c) corresponds to an in- 
crease in entanglement with respect to 0.9 which can be 
as large as 3% for ri = 5.65 and T2 — 37r/2). How- 
ever, as seen in panel (b), the changes in the behaviour 
of the transferred entanglement are much less abrupt (no 
qualitative changes with respect to the passive case of 
Ti = are found up to r2 = 1) and the values of the 
maximally transferred entanglement are lower than the 
corresponding ones obtained for a = 0. The optimality 
of the projection on the excited states of the ancillary 
qubits obtained for a = can be understood at the level 
of the density matrix structure. Indeed, such a projec- 
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tion always results in the postselection of states of qubits 
1 and 2 of the following form: 



gcd,12 = 



f Ai -B\ 

^2 



\~B AiJ 



(9) 



with the density matrix's entries being, in general, very 
complicated functions of ti^2 and r (no qualitative infor- 
mation can be gathered, in general, from the presentation 
of their analytical form). Now, for optimal choices of ti 
an T2, corresponding to the values discussed in Figs.El(c) 
andEI one has Ai,A4,B ^ ^2,3 > so that Qcd,i2 is a 
highly pure state with a very large projection on the Bell 
state (|00) - |ll))/\/2. On the other hand, as soon as 
I a I > 0, contributions different from Eq. ^ appear, be- 
coming increasingly relevant for larger values of jal, thus 
reducing the efficiency of the entanglement transfer. 

The efficiency of the postselected protocol depends on 
the initial squeezing of the CV state as well. Indeed, if we 
consider a small degree of initial entanglement between 
the field modes of the CV system, it can be analytically 
seen that the entanglement between qubits 1 and 2 in the 
postselected protocol described so far is not larger than 
the entanglement transferred through the passive proto- 
col of Refs. 0, SEs'l- In particular, by considering the 
expansion of Eq. |^ into terms up to O(r^) and the joint 
qubits-CV modes dynamics, the postselected density ma- 
trix (a = 1 being in this instance the most favorable case) 
turns out to be equivalent to the non-postselected one. 
Choices of a different from 1 correspond, in this case, to 
a smaller overlap between the reduced state of the qubits 
and (|00) - |ll))/\/2. 

Let us now note that, as a matter of principle, the en- 
tanglement transfer to a pair of two-level systems can be 
further (probabilistically) improved by considering the 
interaction of the entanglement distributor with addi- 
tional pairs of ancillary qubits. To be more specific, 
we have explicitly considered the case of a third pair 
of qubits, labelled with 5 and 6, which interact for a 
rescaled time T3 with modes ai and 02, after the mea- 
surement of the state of qubits 3 and 4 has been accom- 
plished. We assume that the third pair of qubits are also 
found to be in \tjja) iV'a)- It is possible to choose T2 
and T3 so that the maximum of the entanglement trans- 
ferred by this active sequential protocol is even larger 
than what is shown in Fig. |Sl(c). Indeed, for ti ~ 5.65, 
T2 = T3 = TT (and r = 0.86) the transferred entanglement 
reaches ~ 0.99, with a difference from the passive 
transferred entanglement enhanced to 0.093 (see Fig.lSJ. 

This result is per se interesting as it permits to con- 
jecture that the sequentialization of the postselected pro- 
tocol would allow for an asymptotically perfect entangle- 
ment extraction. Moreover, quite remarkably, such an it- 
erative process is efficient, in that a relatively small num- 
ber of iterations does provide an almost perfect trans- 
fer (as we have shown, the results are nearly perfect 





FIG. 7: Entanglement against the rescaled interaction time 
Ti — T for the case |q|^ — 0.95^ and an initial squeezing 
parameter r = 0.86. From (a) to (d) we show the behavior of 
the entanglement function for T2 — 0, T2 — it, T2 — 3tt/2 and 
T2 — 2n. All quantities plotted are dimensionless. 



En 

0.093 \ 



0.0465 



1 2 3 4 5 6 

FIG. 8: Positive part of the difference between the logarithmic 
negativity for n = r and T2 — rs ~ 3 and the passive peak 
value (for which Ej^ ~ 0.9 against the interaction time r. 
An enhancement of about 0.093 is found, which improves the 
result for a single iteration protocol shown in Fig.|S| We have 
considered again r = 0.86 and \af = 0. AH quantities plotted 
are dimensionless. 



already after two iterations). Clearly, with respect to 
the first iteration of the postselection process, the second 
one is much less efficient, in that the entanglement en- 
hancement is a much smaller fraction of the total trans- 
ferred entanglement. This reduction in efficiency is a 
feature characteristic of any 'distillation-like' (in a broad 
sense) process '21^ and sets a practical trade-off between 
the amount of entanglement one would like to extract 
from the distributor and the number of iterations allowed 
by the decoherence time which characterizes the specific 
setup adopted. In many practical cases, like the one we 
are addressing, just two repetitions of the protocol could 
be enough to extract a nearly entire ebit. 
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IV. CONCLUSIONS AND OUTLOOK 

We have considered the transfer and extraction of CV 
entanglement by Jaynes-Cumming interactions with fi- 
nite dimensional atomic-like systems. Our results show 
that microscopic systems or, from a slightly different 
standpoint, systems with Hilbert spaces of 'small' di- 
mensions (< 20) allow for an almost complete extraction 
of the CV entanglement originally stored in a infinite- 
dimensional Hilbert space. On the other hand, the rea- 
sons for the impossibility of a perfect extraction have 
been clearly shown to lie in the feature of the considered 
dynamics, which cannot perfectly mimic that of a quan- 
tum harmonic oscillator. In this respect, the macroscopic 
highly polarised limit (achievable in atomic clouds) can 
be already considered as optimal. However, it could be 
interesting to envisage a realistic microscopic system in 
which the couplings between different levels do actually 
reproduce the harmonic oscillator's ones, described by 
truncated bosonic ladder operators. Such a development 
could deserve further investigation. 



Furthermore, we have shown that a simple and realistic 
sequential postselection protocol would allow for the ar- 
bitrarily perfect extraction of the CV entanglement into 
a Bell state and, even more significantly, that one iter- 
ation of such a protocol would provide with an almost 
perfect symmetric Bell state. Such a strategy could be a 
promising candidate to build up entangled pairs at a dis- 
tance, with remarkable improvements over passive purely 
dynamical strategies. 
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